Calculation of isotope shifts for cesium and francium 
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We perform ab initio calculations of isotope shifts for isotopes of cesium (from A=123 to yl=137) 
and francium (from ^=207 to ^=228). These calculations start from the relativistic Hartree-Fock 
method and make use of several techniques to include correlations. The field (volume) isotope 
shift is calculated by means of an all-order correlation potential method and within the singles- 
doubles partial triples linearized coupled-cluster approach. Many-body perturbation theory in two 
different formulations is used to calculate the specific mass shift. We discuss the strong points 
and shortcomings of the different approaches and implications for parity nonconservation in atoms. 
Changes in nuclear charge radii are found by comparing the present calculations with experiment. 
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I. INTRODUCTION 



Accurate calculations of isotope shifts for many- 
electron atoms are needed to address a number of press- 
ing problems. Although possible changes of isotope abun- 
dances in the early universe are important systematic ef- 
fects in the study of variation of fundamental constants 
0,1^1^, isotope shift data are unavailable for most of the 
spectral lines of interest. Furthermore, studies of parity 
nonconservation (PNC) in atoms have come to the point 
where more information on nuclear structure is needed. 
The analysis of the most precise measurement of PNC 
in cesium reveals that at the present level of accuracy 
(0.35% for experiment, 0.5% for theory, see, e.g. 0) 
there is perfect agreement with the standard model of the 
electroweak interactions. To search for new physics be- 
yond the standard model in low-energy physics, one needs 
to improve the accuracy of the analysis. While improving 
the accuracy of measurements is probably feasible, im- 
proving the accuracy of calculations is problematic. To 
avoid problems with the accuracy of calculations, Dzuba 
et al. [1| suggested that PNC be measured in a chain of 
isotopes. The electronic structure factor would be can- 
celed in the ratio for different isotopes. However, at the 
level of accuracy needed to search for new physics, the 
cancellation is not sufficiently complete. One must ac- 
count for the change in electronic structure caused by 
the change in nuclear charge radius. The change in neu- 
tron distribution presents a separate problem |6|| that will 
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not be discussed further in this work, since the isotope 
shift is not sensitive to the changes in the neutron dis- 
tribution. The change in the nuclear charge radius can 
be obtained by comparing experimental and theoretical 
values of isotope shifts. Moreover, calculations of the iso- 
tope shift provide tests of atomic wave functions at short 
distances, thereby providing another way of evaluating 
the accuracy of atomic PNC calculations. 

We have chosen cesium and francium for our analysis 
since both atoms are being used or planed to be used 
in atomic PNC measurements. Also, both atoms have 
relatively simple electronic structure with one external 
electron above closed shells, making it easier to study 
different theoretical approaches. Accurate experimental 
values of isotope shifts for a range of isotopes of both 
atoms are available. 

Isotope shifts of cesium, francium, and thallium have 
been studied theoretically by A.-M. Martensson-Pendrill 
et al. 0, H, In this paper, we present several dif- 
ferent approaches aimed at improving the accuracy of 
isotope shift calculations and the corresponding analysis 
of changes in nuclear charge radii. We demonstrate that, 
owing to the extremely poor convergence of many-body 
perturbation theory, all-order techniques are needed to 
obtain reliable results. 

In heavy atoms, the isotope shift is dominated by the 
field (or volume) shift (FS), which is the shift in energy 
caused by the change in the nuclear charge radius. In 
light atoms, the isotope shift is dominated by the mass 
shift, which is the change in energy due to differences in 
isotopic masses. Both of these effects are important for 
cesium. The field shift strongly dominates for francium. 
The FS is easier to calculate owing to the simpler form 
of the corresponding operator. We present two different 
all-order methods, both of which use the finite-field ap- 
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proach to reduce the calculation of the FS to a calculation 
of energy. The first method is the all-order correlation 
potential method and the second is the singles-doubles, 
partial triples, linearized coupled-cluster method. In or- 
der to reduce the calculation of the FS to a calculation 
of an energy, we incorporate the rescaled operator of the 
derivative of the nuclear potential with respect to the nu- 
clear radius into the Hamiltonian at each stage of the cal- 
culation. The isotope shift is calculated as the derivative 
of the energy with respect to the scaling parameter. Both 
methods give very close results. Calculations of energies 
and hyperfine structures (hfs) are used as additional tests 
of the accuracy of the calculations. We believe that the 
uncertainty in our calculation of the field shift does not 
exceed 1%. 

For reasons explained below, the same methods cannot 
be used at the present time for the calculation of the mass 
shift. Therefore, we use two different approaches there. 
The first approach is a third-order perturbation theory 
calculation (with certain classes of terms included to all 
orders). The second method which we use in this work 
is a combination of perturbation theory and the finite- 
field approach in which only core polarization diagrams 
are included to all orders. While such a treatment of the 
specific mass shift allowed us to obtain reasonably good 
results suitable for the present work, we stress the need to 
develop a complete all-order technique for the calculation 
of the mass shift, similar to that used here for the field 
shift. 

We extract values of the change in nuclear charge ra- 
dius for different isotopes of cesium and francium by com- 
paring our calculations with available experimental data 
and we discuss the implications of this work for PNC in 
atoms. 



II. METHOD OF CALCULATIONS 

The shift of transition frequency of isotope A' com- 
pared to isotope A can be written as 
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where A and A' are mass numbers of the two isotopes and 
(r^) is the mean square nuclear radius. The first term on 
the right represents the mass shift and the second term 
represents the field shift. The mass shift consists of two 
parts: a normal mass shift (NMS) and a specific mass 
shift (SMS). The normal mass shift constant is expressed 
in terms of the experimental frequency by 



-K'nms = 



1822.888 



(2) 
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56.055 


8 


-1407.65 


34.003 


54.241 


9 


-1404.21 


35.013 


55.172 


10 


-1405.98 


34.494 


54.694 



A. Field shift 

We start with field shift calculations because they are 
easier to carry out and because they play more important 
roles in our analysis owing to their implication for PNC 
in atoms and their dominance for the francium isotope 
shift. 

We use the following form for the FS operator: 



<5V^nuc(r) 



dVnucjr, Rn) 
(IRn 



S{Rn), 



(3) 



where i?Ar is the nuclear charge radius. We assume a 
Fermi distribution for the nuclear charge. The derivative 
dVnuc{r, Rn) /dRN is calculated numerically. 

The change in the nuclear potential ^Vnuc(^) is a small 
perturbation which suggests that perturbation theory 
is probably a proper tool for calculations. However, 
the convergence of perturbation theory in the residual 
Coulomb interaction is extremely poor. In Table ^ we 
present the dependence of the FS constants for the 6s 
and 6p states of Cs on the number of iterations (order of 
perturbation theory) of the linearized Hartree-Fock equa- 
tions for an atom in external field. Here external field 
is produced by the change of nuclear radius and corre- 
sponding equations are equivalent to the random phase 
approximation (RPA) (see, e.g. fTo|). The convergence 
for the 6s state is seen to be very poor while perturba- 
tion theory completely fails for the 6p states, inasmuch as 
one needs more than ten iterations to get stable results. 

Therefore, instead of using perturbation theory, we use 
an all-order finite-field approach similar to our early work 
[Tl|. Calculations of the FS are done for the reference 
isotope A with nuclear potential Vnuc{r, R]sf{A)) replaced 
by 



Vnuc{r,RN{A))+X6Vnuc{r), 



(4) 



while SMS and FS constants -ft^sMS and F are the sub- 
jects of the present calculations. 



where A is scaling parameter. The value of A is chosen in 
such a way that the corresponding change in the nuclear 
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TABLE II: Field shift constants F (MHz/fm^) for Cs and Fr 
in different approximations. 

Cs Fr 





6s 


6pi/2 


6P3/2 


7pi/2 


7s 


7pi/2 


7p3/2 


HF 


-1270 


-15.7 


0.0 


-5.6 


-14111 


-458 


0.0 


RPA 


-1405 


34.6 


54.9 


12.4 


-15819 


-209 


510 




-2050 


17.4 


45.4 


8.7 


-22358 


-697 


313 


B0(E(2>) 


-2119 


17.6 


46.7 


10.0 


-22447 


-759 


301 


BO(E°°) 


-1914 


22.4 


51.2 


10.6 


-20463 


-693 


303 


SD+E3 


-1894 


23.9 


52.2 


11.3 


-20188 


-640 


361 


Otlier 


-2069" 


42.78" 


70.53" 


15.17" 


-20782' 


-696' 


245* 



"Hartley and Martensson-Pendrill 
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potential is much larger than the numerical uncertainty 
of the calculations but is still sufficiently small for the 
final energy to be a linear function of A. The PS constant 
for a particular atomic state v is then found as 

This approach reduces the calculation of the FS to a cal- 
culation of energy. We use two different techniques to 
calculate energies. One is the all-order correlation po- 
tential method 0, ^| (also called perturbation theory 
in the screened Coulomb interaction). Another is the 
singles-doubles linearized coupled cluster method com- 
bined with many-body perturbation theory (MBPT) to 
account for missing third-order diagrams (SD-I-E3). 

We present the results in Table UTI in order of improv- 
ing approximations. The first line (HF) gives the average 
value of the 6Vnuc{r) over Hartree-Fock wave functions. 
The second line (RPA) gives the result of Hartree-Fock 
iterations with the potential given by Eq. Q. The re- 
sult of this calculation (linear in the scaling parameter 
A) is equivalent to the RPA. The next line ((S*^^^)) in- 
cludes second-order correlation corrections by means of 
many-body perturbation theory. Since we use "dressed" 
basis states (states calculated in a potential given by the 
Eq. with the FS operator included) to calculate these 
corrections, all third-order terms are included (first in 
FS and second in Coulomb interaction) as well as chains 
of higher-order terms corresponding to core polarization. 
We use the notation E for the correlation correction op- 
erator (correlation potential) . Details of the use of S in 
atomic calculations can be found elsewhere ^3 • The line 
labeled BO(S^^^) presents results obtained by including 
the operator E^^^ in the Hartree-Fock equations for the 
valence electron and calculating Brueckner orbitals (BO) 
and the corresponding energies. These results differ from 
those in the previous line by higher-order contributions 
in t (E2,E^ etc.). Finally, the line BO(E°°) presents 
results obtained with the all-order E, in which screening 
of Coulomb interaction and hole-particle interactions are 



included in E°° to all orders (see, e.g. for details) and 
E°° is used to calculate BO's. These are the most accu- 
rate results obtained in the all-order correlation potential 
(CP) method. We compare them with results obtained 
in the SD-I-E3 approach presented in next line. We re- 
fer the reader to Ref. jl^l for the details of the (SD-I-E3) 
all-order energy calculation. For the 6s — 6p and 6s — 7p 
intervals (the only ones important for the analysis of the 
experimental results), the two methods agree to about 
1%. 

In the last line of Table ^1 we present the FS calcu- 
lations of Hartley and Martensson-Pendrill 0, |1] . They 
use an approximation very similar to that listed on line 
(E^^)). Naturally, the results are also very close. The 
larger discrepancy seen for the 6p states of Cs is proba- 
bly due to larger contributions from structural radiation 
diagrams (which are called internal diagrams in Ref. 0). 
In our approach, calculation of the IS is reduced to a cal- 
culation of energy and no special treatment for the struc- 
tural radiation diagrams is needed. In contrast. Hartley 
and Martensson-Pendrill evaluated these diagrams using 
a procedure based on a modification of the basis. We be- 
lieve that our calculations are more accurate since (a) the 
finite-field approach ensures that no important diagrams 
are missed, (b) we have very good agreement between 
two very different methods, and (c) we have very good 
agreement with experiment for both energies and hyper- 
fine structures (see below). 

The final result of Hartley and Martensson-Pendrill 
for the 6s state of Cs (-2000 MHz/fm^) was obtained 
by rescaling the ab initio value using a comparison be- 
tween theoretical and experimental hyperfine structure 
(hfs) constants. Their calculated value for the hfs of 6s 
is larger than the experimental value by 3%. Therefore, 
they reduced the FS constant by the same ratio. Since 
we include higher-order correlations which bring calcu- 
lated values of the hfs constants and the energies into 
very good agreement with experiment, we can check how 
well this rescaling works. 

Table lHll presents the data for the second- and all-order 
energies, hfs constants, and FS constants for Cs and Fr. 
We see that higher-order correlations reduce the ground 
state energy by 3%, the hfs constant by 5 to 8% and FS by 
9 to 10%. Therefore, energies cannot be used for rescal- 
ing of the FS constants, and accuracy of extracting the 
higher-order correlation correction to the FS by rescaling 
the second-order correction using hfs data is between 20 
and 40%. This corresponds to 2 to 4% accuracy in final 
results for s-states. Moreover, no reliable rescaling can 
be done for p-states. 

Note that the rescaling of FS is very different from 
the "fitting of the energy" procedure used in some of 
our works (see, e.g. |j|) and in the next section. In 
that procedure, we replace E in equations for valence 
BO's by the rescaled operator AE with A chosen to fit 
experimental energies exactly. New BO's are then used 
to calculate matrix elements. This procedure changes 
energies and matrix elements at different rates usually 
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TABLE III: Contributions of the higher-order correlations 
to the energies, hyperfine constants (hfs), and field shift con- 
stants (FS) of Cs and Fr. 



BO(St^>) BO(S°°) Ratio Exp 

Cesium energies (cm~^) 



6s 


32375 


31470 0.97 


31407" 


6pi/2 


20524 


20296 0.99 


20229" 


6p3/2 


19926 


19728 0.99 


19675" 




133 


Cs hfs (p/=0.7377208) (MHz) 


6s 


2459 


2270 0.92 


OOOQ 0^ 

zzyo.z 


6pi/2 


314 


295 0.94 


291.89' 


6P3/2 


51.8 


48.7 0.94 


50.275"^ 






Cesium FS (MHz/fm^) 




6s 


-2119 


-1914 0.90 




6Pl/2 


17.6 


22.4 1.3 




6p3/2 


46.7 


51.2 1.1 








Francium energies (cm~^) 




7s 


34089 


32899 0.97 


32849' 


7Pl/2 


20986 


20711 0.99 


20612' 


7p3/2 


19164 


18976 0.99 


18925' 




2"Fr hfs (p/=0.888) (MHz 


) 


7s 


9269 


8769 0.95 


8713.9(8/ 


7pi/2 


1261 


1193 0.95 


1142.0(3)« 


7P3/2 


98.3 


102.4 1.04 


94.9(3)^ 






Francium FS (MHz/fm^) 




7s 


-22447 


-20463 0.91 




7Pl/2 


-759 


-693 0.91 




7P3/2 


301 


303 1.01 





"Moore [H 

''Arimondo et al. Ilfl 

'^Rafac and Tanner lldl 

■^Tanner and Weiman llTl 

'Bauche et al. ^3, Arnold et al. |1H|I 

^Ekstrom et al. 12(1 

^Grossman et al. |2lj 

bringing the later into significantly better agreement with 
experiment. 

The comparison of final energies and hfs constants with 
experiments presented in Table ITTll is a further indication 
that the accuracy of the present calculations of FS con- 
stants is at the level of 1%. 



B. Specific Mass Shift 

The finite-field approach used in the previous section 
to calculate the field shift can also be used to calculate 
the specific mass shift. To do so we need to redefine the 
Coulomb interaction in the following way 

{ah\e^/ri2\cd) (a5|eVri2|cd) + X{ab\pi ■ P2\cd), (6) 

where A is a scaling parameter and p is the electron mo- 
mentum (see ^2 for details). 



The substitution © can be easily done in the Hartree- 
Fock approximation or in perturbation theory calcula- 
tions. However, the methods used in the previous section 
to include higher-order correlation corrections in the FS 
constant are not applicable here. Note that for FS calcu- 
lations only the Hartree-Fock program needs to be mod- 
ified to incorporate the change in nuclear potential 
By contrast, for SMS calculations ewery program must be 
modified. While this is still straightforward in Hartree- 
Fock and perturbation theory codes, it becomes much 
more difficult in higher orders. In the correlation poten- 
tial method, inclusion of higher-order correlations is done 
by summing the matrix geometric progression jl^ 

Q = Q + QnQ + QIlQIlQ + ... (7) 

where Q is the Coulomb interaction ((9fc(ri,r2) = 
r^/r>^^), n is polarization operator, and Q is screened 
Coulomb interaction (see, e.q. [13 for details). To in- 
clude the SMS operator in this summation, one would 
need to modify Coulomb interaction in the following way 

Q^Q + XP, (8) 

where P is the coordinate representation of the SMS op- 
erator, which would lead to the correct expressions for 
radial integrals © when integrated over the wave func- 
tions. It is clear that there is no such representation 
for the SMS operator. This problem does not appear in 
the SD-I-E3 method, since everything there is expressed 
in terms of Coulomb integrals which can be modified ac- 
cording to Eq. © . While such modification is as straight- 
forward as for perturbation theory codes, technically it is 
not an easy task. The problem here is not only with the 
large number of terms which must be modified, but also 
with different symmetry properties of the SMS operator. 
Exchanging indexes a and c (or b and d) in Eq. © leaves 
the Coulomb part of the equation unchanged while the 
SMS contribution changes sign! While we stress that it 
would be extremely useful to have a finite-field program 
for the SMS, we must leave this for future work. 

In the present work, we use two less sophisticated (and 
less accurate) approaches. The first is perturbation the- 
ory and the second is a finite-field approach in which 
only the second-order correlation operator E^^^ is used. 
In the perturbation theory calculation, we express the 
SMS operator P = J2i<j Pi ' Pj ^ sum of a normally 
ordered one-particle operator S and normally ordered 
two-particle operator T. We carried out an all-order cal- 
culation of the matrix element of S; the calculation of 
the matrix element of T is complete through third-order 
(first order in the SMS operator and up to second order 
in Coulomb interaction) . The results of the perturbation 
theory calculation for the SMS for Cs and Fr are pre- 
sented in Table HVl The lowest-order values are given in 
row S**-^-*, the results of the all-order singles-doubles cal- 
culation of the matrix element of S are given in row SD, 
and the second- and third-order matrix elements of T are 
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TABLE IV: Perturbation theory contributions to specific 
mass isotope shift constants (GHZ amu) in Cs and Fr 







Cs 






Fr 






6s 


GPl/2 


6P3/2 


7s 


7Pl/2 


7P3/2 




-781.3 


-191.4 


-168.7 


-1359.7 


-260.0 


-187.0 


SD 


316.5 


139.9 


140.3 


363.7 


182.7 


182.4 




286.8 


58.5 


51.2 


499.9 


78.9 


57.5 


rp(3) 


-136.2 


-30.1 


-27.8 


-296.0 


-53.7 


-45.0 


Total 


-314.2 


-23.1 


-5.0 


-786.1 


-53.0 


7.9 



given m rows labeled T^^) ^nd T^^). The total values of 
the specific mass shift constants are listed in the last row. 

Table presents results of the finite-field approach. 
The line HF gives the expectation values of the SMS op- 
erator over HF wave function of the valence electron, ft 
is equivalent to the line S'^^^ of Table HVl The difference 
in numerical values is due to the fact that the relativistic 
form of the momentum operator was used in the PT cal- 
culations (S**-^^) while the non-relativistic operator was 
used in the finite-field calculation. Note that relativistic 
corrections for s ~ p intervals do not exceed 4% for Cs 
and 8% for Fr. 

This is a negligible contribution since the mass shift 
itself is small for heavy atoms owing to the huge sup- 
pression by the mass factor (see Eq. (QJ). Note, however, 
that relativistic corrections are probably very important 
for highly-charged ions. 

The line RPA in Table IVI presents results of HF itera- 
tions with the SMS operator included in the HF poten- 
tial by redefining of the Coulomb interaction according 
to Eq. (O. The line S^^^ includes second-order correla- 
tion corrections. Finally, the line Brueck presents results 
for valence electron Brueckner orbitals calculated with 
using We have also included two ways of simulat- 

ing higher-order correlations to try to further improve 
the results. Note that Brueckner orbitals with second 
order S considerably overestimate the correlation correc- 
tion to the energy. They probably have the same effect 
on the SMS. Therefore, we reduce the total correlation 
correction to the SMS in two different ways. 

Firstly, we note that iterations of E^^-* enhance cor- 
relation corrections to the SMS by a larger factor than 
for energies. If we use the energy ratio to determine the 
enhancement instead, the total correlation correction to 
the SMS is smaller. The corresponding interpolated re- 
sults are presented in the line Interp in Table E] Note 
that the correction is huge. It can even change sign of 
the SMS constant. This is very different from the FS con- 
stants discussed in previous section. For the FS constants 
higher-order correction is small and should be treated ac- 
curately. This is why rescaling works for s-states only. 

Secondly, we rescaled the operator S(2) while calculat- 
ing Brueckner orbitals to fit the experimental energies. 
Scaling factors are A(6s) = 0.802, A(6p) = 0.85 for Cs 
and A(7s) = 0.786, \{7p) = 0.85 for Fr. This procedure 



also reduces the correlation correction to the SMS. Cor- 
responding results are presented in line Fit of Table Ivl It 
is interesting that the two procedures give close results. 

Comparison of results in Table IIVI and Table re- 
veals the poor convergence of perturbation theory and 
the significant difference in final results between the two 
calculation methods. Note that the two methods are 
equivalent at second order in the Coulomb interaction. 
The difference comes from higher orders. Perturbation 
theory calculations use the SMS matrix elements which 
are just expectation values of the SMS operator over HF 
wave functions (HF matrix elements). By contrast, the 
finite-field approach corresponds to including "dressed" 
SMS matrix elements in which certain chains of Coulomb 
diagrams are included to all orders by iterating the HF 
equations (RPA matrix elements). By comparing lines 
HF and RPA of TablelVlone can see that this indeed must 
lead to large differences in final results. Note that the 
Brueckner-orbital calculations are also in better agree- 
ment with the calculations of Hartley and Martensson- 
Pendrill HQ. 

To check whether Brueckner orbitals really give better 
results than PT for the SMS, we have performed calcu- 
lations for potassium and rubidium. For these atoms, 
"experimental" values of the SMS can be inferred by 
subtracting the NMS and FS from known experimen- 
tal values of the isotope shifts. The NMS is given by 
Eq. |(2Jl while the field shift is calculated as F5{r'^). The 
field shift constant F can be calculated to high precision 
as described in the previous section. The corresponding 
values are presented in Table IVll Values of 5{r'^) be- 
tween the most abundant isotopes are also known from 
experimental studies Table IVTll presents the extrac- 
tion of the "experimental" SMS between "^^K and "^^K 
for the 4s — Api/2 transition and between ®^Rb and *^Rb 
for the 5s — 5p3/2 transition. In Table IVlTTl we com- 
pare these experimental results with calculations. One 
can see that interpolated, fitted, and unfitted Brueckner 
orbital results are closest to experiment while the pertur- 
bation theory results are much farther away. It is natural, 
therefore, to take the average value of the most accurate 
results {Brueck, Interp and Fit) as a final central point 
of the calculations, while using differences between these 
results as an estimate of numerical uncertainty. 



III. RESULTS AND DISCUSSION 

Final values for the mass and field shift constants are 
presented in Table |IX| These are the BO (S°°) results 
for FS (see Table HH) and the average of Brueck, Interp 
and Fit results for SMS (see Table W\ . We use these 
values to analyze experimental data and to extract the 
change in nuclear charge radius for a range of isotopes of 
cesium and francium. The results are presented in Tables 
IXland lXll We present two uncertainties for the (5(r^) for 
Cs. The first one is experimental and the second one is 
theoretical. The theoretical uncertainty is dominated by 
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TABLE V: SMS constants for Cs and Fr in different approximations (GHZ amu). 





6s 


6pi/2 


Cs 

6P3/2 


7pi/2 


7P3/2 


7s 


Fr 

7pi/2 


7p3/2 


HF 


-773.18 


-208.70 


-170.40 


-74.68 


-61.44 


-1330.48 


-317.79 


-190.87 


RPA 


-355.55 


-76.96 


-40.31 


-27.59 


-14.56 


-666.71 


-127.90 


-10.48 


S(2) 


-133.19 


-22.83 


20.66 


-5.63 


8.81 


-334.15 


-85.24 


60.08 


Brueck 


6.75 


6.85 


47.33 


0.95 


14.31 


-110.29 


-45.11 


89.86 


Interp 


-99.84 


-14.39 


29.38 


-4.31 


10.21 


-288.26 


-77.90 


69.96 


Fit 


-89.31 


-10.26 


30.17 


-2.94 


10.47 


-311.96 


-63.83 


70.35 


prpa 


-314.2 


-23.1 


-5.0 






-786.1 


-53.0 


7.9 


Other 


-23.5'' 


-36.6' 


9.2'' 






-570' 


-154' 


-18' 



"From Table ITVI 

'Hartley and Martensson-Pendrill 0| 
'^Martensson-Pendrill Q 



TABLE VL Field shift constants for K and Rb (MHz/fm^) 

K 4s K 4pi/2 Rb 5s Rb 5p3/2 

-104.20 4.04 -551.85 15.60 



TABLE VIL Extracting an "experimental" SMS for K and 
Rb. 



IS(exp) 


5{r^r 


FS 


NMS 


SMS 


(MHz) 


(fm^) 


(MHz) 


(MHz) 


(MHz) 




39- 


-«K 4s- 


4pi/2 




235.25(75)* 


0.105 


-11.33 


267.11 


-20.53(75) 




85- 


^^■'Rb 5s - 


5p3/2 




77.992(20)' 


-0.042 


23.83 


57.00 


-2.84(2) 


"Angeli I2I 



''Touchard et al. El 
'^Banerjee et al. |2fi| 

the uncertainty in the SMS constant. Strong cancellation 
between the field shift and the normal mass shift makes 
these results sensitive to the SMS. However, the very 
poor convergence of perturbation theory make it difhcult 
to predict the SMS to high accuracy. We stress once more 
the need for accurate all-order techniques to calculate the 
SMS. 

It has been long known (see, e.g. |2^) that the change 
of nuclear radius along the chain of Cs isotopes is slower 
than expected from the formula Rn — 1.1^^/"^. One 
possible explanation for this fact is that neutrons and 
protons have different distributions. According to data 
deduced from antiprotonic atoms [23|, the radius of the 
neutron distribution for ^'^^Cs is by 0.16 ± 0.06 fm larger 
than radius of the proton distribution. It is interest- 
ing to note that under certain assumptions, the isotope 
shift data is in very close agreement with the antiprotonic 
atom data. These assumptions are: (a) the A^^^ law is 
still valid but for the total nuclear radius (including the 



TABLE VIII: Finite-field SMS constants for K and Rb in 



different approximations (GHZ amu). 







K 






Rb 






4s 


4pi/2 


A 


5s 


5P3/2 


A 


HF 


-387.27 


-120.56 


-266.71 


-587.21 


-144.92 


-442.29 


RPA 


-193.62 


-61.51 


-132.11 


-263.55 


-45.23 


-218.32 




-51.71 


-25.95 


-25.76 


-69.46 


2.45 


-71.91 


Brueck 


-4.36 


-15.00 


10.64 


14.73 


19.25 


-4.52 


Interp 


-37.52 


-22.39 


-15.13 


-46.17 


7.69 


-53.86 


Fit 


-29.08 


-19.07 


-10.01 


-31.81 


11.85 


-43.66 




-74.7 


-24.7 


-50.0 


-163.3 


-11.8 


-151.5 


Exp 






-16.4 






-10.5 



"Safronova and Johnson |22| 



TABLE IX: Final values of the IS shift constants for Cs and 
Fr used for the analysis of the experimental data 





Transition 


KsMS 


-R'nms 


F 






GHz amu 


GHz amu 


MHz/fm^ 


Cs 


6s - 6p3/2 


-96(56) 


192.9 


-1965(20) 


Cs 


6s - 7pi/2 


-59(65) 


357.9 


-1925(20) 


Fr 


7s - 7p3/2 


-314(113) 


229.0 


-20766(208) 



"neutron skin") not just charge radius; (b) the neutron 
and proton distributions are very close for the neutron 
poor nuclei. This is one more argument in favor of the 
"neutron skin" correction to the PNC in Cs "2^. 

Table IXII presents IS data for Fr. The experimental 
accuracy for the IS is very high. Therefore, only the the- 
oretical error is presented for the 5{r'^). The SMS is less 
than 1% of the total IS for Fr. Also, there is a strong can- 
cellation between the SMS and the NMS, which makes 
the extracted values of the change in nuclear charge ra- 
dius insensitive to the SMS. The theoretical uncertainty 
is therefore quite low and at the level of 1% which is 
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TABLE X: Isotope shift (MHz) and change of nuclear radius 
(fm^) between ^^^Cs and other cesium isotopes 



TABLE XL Isotope shift (MHz) and change of nuclear radius 
(fm^) between ^^^Fr and other francium isotopes 



A 6s- 


SMS 


NMS 


FS 


Exp IS" 




A 


SMS 


NMS 


FS 


Exp IS" 


y j 


123 7pi/2 

^ ' IT i- j ^ 


35.88 


-218.82 


441.94 


259(12) 


-0.230(10)(36) 


207 


35.73 


-26.09 


5229.36 


5239(4) 


-0.252(3) 


124 Tdi /o 


32.03 


-195.35 


424.31 


261(6) 


-0.220(5) (30) 


208 


28.44 


-20.77 


4995.33 


5003(3) 


-0.241(2) 


125 7pi/2 


28.25 


-172.25 


296.01 


152(11) 


-0.154(11)(32) 


209 


21.23 


-15.50 


3127.27 


3133(2) 


-0.151(2) 


126 7pi/2 


24.52 


-149.53 


333.01 


208(7) 


-0.173(6)(23) 


210 


14.09 


-10.29 


2599.20 


2603(1) 


-0.125(1) 


127 7pi/2 


20.85 


-127.16 


200.30 


94(13) 


-0.104(14)(26) 


211 


7.01 


-5.12 


899.11 


901(3) 


-0.0439(4) 


128 7pi/2 


17.24 


-105.14 


242.90 


155(6) 


-0.126(5)(16) 


213 


-6.94 


5.07 


-1639.13 


-1641(2) 


0.0790(8) 


129 7pi/2 


13.69 


-83.46 


122.77 


53(9) 


-0.064(11)(18) 


220 


-53.78 


39.28 


-20792.29 


-20806.8(0.5) 


1.001(10) 


130 7pi/2 


10.19 


-62.11 


107.93 


56(8) 


-0.056(8)(11) 


221 


-60.23 


43.99 


-23553.76 


-23570(2) 


1.134(11) 


130 6p3/2 


16.73 


-33.48 


90.75 


74.0(2.2) 


-0.046(1)(5) 


222 


-66.62 


48.66 


-26244.03 


-26262(3) 


1.264(13) 


131 7pi/2 


6.74 


-41.09 


25.35 


-9(6) 


-0.013(8)(11) 


223 


-72.96 


53.28 


-27902.32 


-27922(2) 


1.344(14) 


131 6p3/2 


11.07 


-22.15 


21.48 


10.4(1.6) 


-0.011(2)(7) 


224 


-79.24 


57.87 


-30869.63 


-30891(1) 


1.487(15) 


132 7pi/2 


3.34 


-20.39 


77.05 


60(15) 


-0.040(10) (3) 


225 


-85.46 


62.41 


-32273.95 


-32297(1) 


1.554(16) 


132 6p3/2 


5.49 


-10.99 


79.80 


74.3(1.3) 


-0.041(7)(2) 


226 


-91.62 


66.91 


-34376.29 


-34401(1) 


1.655(17) 


134 6p3/2 


-5.41 


10.83 


27.68 


33.1(2.5) 


-0.014(1)(1) 


227 


-97.74 


71.38 


-38325.64 


-38352(2) 


1.846(19) 


135 6p3/2 


-10.74 


21.49 


-47.15 


-36.4(2.0) 


0.024(1)(2) 


228 


-103.79 


75.80 


-40049.01 


-40077(5) 


1.929(20) 


137 7pi/2 
137 6p3/2 


-12.89 
-21.17 


78.58 
42.36 


-169.70 
-168.59 


-104(6) 
-147.4(2.5) 


0.088(5)(12) 
0.086(1)(7) 


"Coc ei al. IsntlsH 









"Huber ei al. Thibault et al. |29| 



mostly the uncertainty in the FS constant. Our results 
for S{r^) differ from those obtained in @ by about 2.4%. 
This is due to the difference in the FS constant in the 
7s — 7p3/2 transition. Our value is -20.8 GHz amu while 
a value of -21.0 GHz amu is used in @. Since 2% un- 
certainty is estimated in Q, and we believe that our ac- 
curacy is 1%, we can say that the results are in good 
agreement. 

Since nuclear radii change quite significantly along the 
chain of francium isotopes, it is important to check how 
PNC matrix elements are affected by this change. We 
have conducted a numerical test for the change of the ma- 
trix element W{Rn) = {7s\H™^\7pi/2) with the change 
of nuclear radius, assuming that proton and neutron dis- 
tributions remain the same. The numerical results can 
be presented in a form 



WjRN) 
W{Ro) 



1 -0.21( 



R 

Rn 



N 



1), 



(9) 



where Rq is the nuclear radius of a reference isotope. 
Note that total change between the lightest and the heav- 
iest isotopes in Table IXll is almost 1%. This is a signif- 
icant change and should be taken into account in any 
future analysis of the PNC in the chain of Fr isotopes. 

Let us now discuss the role that isotope shift calcu- 
lations may play in the study of PNC in atoms. As it 
is clear from the results presented above, for atoms in 
the middle of the periodic table (Cs, Ba, etc.) the in- 
formation that can be extracted from the IS calculations 
and used for the PNC analysis is limited by the accuracy 
of the SMS calculations. Unless adequate methods are 



developed to significantly improve the accuracy of such 
calculations it is unlikely that the IS data will provide 
information of any practical use for PNC analysis. 

The situation is very different for heavier atoms, such 
as Fr, Tl, etc. The mass shift is small and the corre- 
sponding uncertainty can be reduced to an acceptable 
level. On the other hand, calculations of the FS are much 
easier and can be done with accuracy of about 1% or pos- 
sibly better. There are several possibilities arising from 
this fact. To use the IS data to test the electron wave 
function one needs to know the value d{r'^) from an in- 
dependent source. Such data can be obtained from elec- 
tron scatterin g, a nalysis of X-rays from muonic atoms, 
etc. (see, e.g. [23). However, the accuracy of that data 
is often insufficient for the PNC purposes. One can do 
a consistency test instead. If the isotope shift is known 
for several different transitions, one can check whether 
comparison of the theory with experiment leads to the 
same value of (5(r^) for all transitions. 

There is another possibility for many-electron atoms. 
If the IS is known for an ion with simple electronic struc- 
ture (one electron above closed shells), then calculations 
of the IS for this ion can be used to extract the value of 
(5(r^). Because of simple electronic structure, the calcu- 
lations are relatively simple and can be done very accu- 
rately. Then this value of (5(r^) can be used in the IS and 
PNC analysis for a neutral atom. 



IV. CONCLUSION 

We have developed methods to calculate the isotope 
shift for many-electron atoms with one external electron 
above closed shells. While methods for the field shift 
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seem to be adequate and capable of producing results at 

the 1% or better level of accuracy, methods for the SMS 
need further consideration. It would be useful to have 
an all-order technique similar to that used in the FS cal- 
culations to address the problem of the very poor con- 
vergence of perturbation theory. We use our calculations 
for cesium and francduni to extract the change in nuclear 
charge radius for chains of isotopes in both atoms. We 
have demonstrated that, at least for heavy atoms, calcu- 
lations and measurements of the isotope shifts may pro- 
vide important information for the analysis of the PNC 



in atoms. 
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